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Olivine-type  LiFeP04  is  widely  considered  as  a  candidate  for  Li-ion  battery  electrodes,  yet  its  applicability 
in  the  pristine  state  is  limited  due  to  poor  ionic  and  electronic  conduction.  Doping  can  be  employed  to 
enhance  the  material’s  electrical  conductivity.  However,  this  should  be  understood  as  incorporating  elec¬ 
trically  active  impurities  to  manipulate  the  concentration  of  native  point  defects  such  as  lithium  vacancies 
and  small  hole  polarons  which  are  responsible  for  ionic  and  electronic  conduction,  respectively,  and  not 
as  generating  band-like  carriers.  Possible  effects  of  monovalent  (Na,  K,  Cu,  and  Ag),  divalent  (Mg  and  Zn), 
trivalent  (Al),  tetravalent  (Zr,  C,  and  Si),  and  pentavalent  (V  and  Nb)  impurities  on  the  ionic  and  elec¬ 
tronic  conductivities  of  LiFeP04  are  analyzed  based  on  results  from  first-principles  density-functional 
theory  calculations.  We  identify  impurities  that  are  effective  (or  ineffective)  at  enhancing  the  concen¬ 
tration  of  lithium  vacancies  or  small  hole  polarons.  Based  on  our  studies,  we  discuss  specific  strategies 
for  enhancing  the  electrical  conductivity  in  LiFeP04  and  provide  suggestions  for  further  experimental 
studies. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Olivine-type  LiFeP04  has  been  proposed  as  a  candidate  for 
rechargeable  Li-ion  battery  electrodes  because  of  its  structural  and 
chemical  stabilities,  high  intercalation  voltage,  high  theoretical  dis¬ 
charge  capacity,  environmental  friendliness,  and  potentially  low 
cost  [  1  -3  ].  However,  poor  ionic  and  electronic  conduction  are  major 
challenges.  In  addition  to  nanostructuring  and  carbon  coating,  dop¬ 
ing  has  been  considered  as  an  important  path  toward  enhancing 
the  electrical  conductivity.  Since  it  was  reported  by  Chung  et  al. 
[4]  that  doping  LiFeP04  with  impurities  such  as  Mg,  Ti,  Zr,  and  Nb 
significantly  enhances  the  conductivity,  there  have  been  numer¬ 
ous  experimental  works  carried  out  along  this  path.  Nonetheless 
the  role  of  dopants  in  the  conductivity  enhancement  is  still  under 
debate  [5,6].  Thus  it  is  important  to  understand  the  mechanisms 
for  ionic  and  electronic  conduction  and  the  effects  of  impurities  on 
the  ionic  and  electronic  conductivities.  Such  an  understanding  is 
essential  in  formulating  strategies  for  improving  the  electrical  (i.e., 
ionic  and/or  electronic)  conductivity  and  hence  the  electrochemical 
performance  of  the  material. 


*  Corresponding  author  at:  Center  for  Computational  Materials  Science,  Naval 
Research  Laboratory,  4555  Overlook  Ave.  SW,  Washington,  DC  20375,  USA. 

Tel.:  +1  202  404  8626. 

E-mail  address:  michelle.johannes@nrl.navy.mil  (M.D.  Johannes). 

0378-7753 /$  -  see  front  matter.  Published  by  Elsevier  B.V. 
doi:  10.1 01 6/j.jpowsour.201 2.01. 126 


Recently,  we  have  carried  out  comprehensive  first-principles 
density-functional  theory  (DFT)  studies  of  native  point  defects  and 
defect  complexes  in  LiFeP04  [7].  Based  on  a  detailed  analysis  of 
the  structure,  energetics,  and  migration  of  the  native  defects,  we 
arrived  at  the  following  main  conclusions:  (i)  native  point  defects 
such  as  small  hole  polarons  (p+),  negatively  charged  lithium  vacan¬ 
cies  {V~ ),  negatively  charged  lithium  antisites  (Lipe),  and  positively 
charged  iron  antisites  (Fe+)  have  low  formation  energies  and  hence 
are  expected  to  be  present  in  the  material.  The  mobility  of  Fe+  is  low 
compared  to  V~,  supporting  earlier  suggestions  that  Fe+  impedes 
Li  diffusion  and  hence  reduces  the  electrochemical  activity  [8,9]. 
The  relative  concentrations  of  these  defects  are,  however,  sensi¬ 
tive  to  the  experimental  conditions  during  synthesis.  This  suggests 
that  one  can  suppress  or  enhance  certain  native  defects  via  tuning 
the  synthesis  conditions,  (ii)  Native  defects  in  LiFeP04  cannot  act  as 
sources  of  band-like  electrons  and  holes,  and  the  material  cannot  be 
doped  n-  orp-type.  Any  attempt  to  deliberately  shift  the  Fermi  level 
to  the  valence-band  maximum  (VBM)  or  conduction-band  mini¬ 
mum  (CBM),  e.g.,  via  doping  with  acceptors  or  donors,  will  result  in 
spontaneous  formation  of  compensating  native  defects  that  coun¬ 
teract  the  effects  of  doping,  (iii)  The  ionic  conduction  occurs  via 
diffusion  of  V~  whereas  the  electronic  conduction  proceeds  via 
hopping  of  p+,  which  confirms  earlier  suggestions  of  a  polaronic 
mechanism  in  the  electronic  conduction  [10-13].  The  ionic  con¬ 
duction  is  effectively  one-dimensional  along  the  b-axis,  assuming 
there  are  no  other  native  defects  or  extrinsic  impurities  with  low 
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VBM  HT  CBM 


Fig.  1.  Formation  energies  of  intrinsic/native  defects  Vu  and  p+  (solid  blue  lines) 
and  of  a  donor-like  impurity  (solid  red  line),  plotted  as  a  function  of  Fermi  level 
with  respect  to  the  VBM.  In  the  absence  of  extrinsic  electrically  active  impurities, 
the  Fermi  level  of  LiFeP04  is  at  /4nt  where  the  formation  energies  and  thus  concen¬ 
trations  of  V~  and  p+  are  equal.  When  the  donor-like  impurity  is  incorporated  into 
the  material  with  a  concentration  higher  than  that  of  the  positively  charged  native 
defects,  the  Fermi  level  is  shifted  toward  the  CBM,  lowering  (increasing)  the  forma¬ 
tion  energy  of  V~  ( p+ ).  As  a  result,  the  ionic  (electronic)  conductivity  will  be  enhanced 
(reduced).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  the  article.) 

mobility  that  block  the  Li  channels;  whereas  the  electronic  conduc¬ 
tion  is  effectively  two-dimensional  in  the  b-c  plane. 

It  is  evident  from  the  conduction  mechanisms  of  LiFeP04  that 
the  ionic  and  electronic  conductivities  are  dependent  on  mobility 
and  concentration,  and  hence  on  the  migration  barriers  and  forma¬ 
tion  energies,  of  V~  and  p+.  This  opens  the  door  to  manipulating 
the  ionic  and  electronic  conductivities  of  the  material  via  manip¬ 
ulating  defect  concentrations.  In  fact,  LiFeP04  samples  with  high 
concentrations  of  V~  and/or  p+  (and  a  low  concentration  of  Fe+) 
can  be  obtained  in  experiment  under  appropriate  synthesis  condi¬ 
tions  [7].  Also,  despite  the  fact  that  LiFeP04  cannot  be  doped  n-  or 
p- type,  one  can  still  incorporate  into  the  material  electrically  active 
impurities  that  are  effective  in  shifting  the  position  of  the  Fermi 
level  away  from  that  determined  by  the  intrinsic  native  defects, 
slightly  toward  the  CBM  (for  impurities  that  exhibit  a  donor-like 
effect)  or  VBM  (for  impurities  that  exhibit  an  acceptor-like  effect). 
This  results  in  lowering  the  formation  energy  and  hence  enhancing 
the  concentration  of  V~  or  p+.  Thus  “doping”  in  LiFeP04  should  be 
understood  as  a  way  to  manipulate  the  concentration  of  the  native 
defects  that  are  responsible  for  ionic  and  electronic  conduction, 
rather  than  as  generating  band-like  carriers. 

The  mechanism  for  manipulating  the  concentrations  of  p+  and 
V~  by  incorporating  a  donor-like  impurity  is  illustrated  in  Fig.  1 . 
Note  that,  in  the  absence  of  extrinsic  electrically  active  impurities, 
the  Fermi  level  of  insulating  materials  such  as  LiFeP04  is  at  /xe  = 
/4nt,  determined  by  a  charge  neutrality  condition  that  involves 
all  possible  intrinsic/native  defects.  This  condition  requires  that 
defects  with  different  charge  states  coexist  in  a  proportion 
that  maintains  overall  charge  neutrality  in  the  material  [14-17]. 
Because  of  the  exponential  dependence  of  the  concentration  on 
defect  formation  energy  (see  the  next  section),  /4nt  is  predomi¬ 
nantly  determined  by  the  positively  and  negatively  charged  defects 
with  the  lowest  formation  energies  (in  the  illustrated  example, 
these  defects  are  V~  and  p+).  If  electrically  active  impurities  are 
incorporated  into  LiFeP04  with  concentrations  higher  than  that 
of  the  native  defects,  the  Fermi  level  may  be  shifted  away  from 


/4nt  as  the  charge  neutrality  condition  is  re-established  [14,15],  as 
illustrated  in  Fig.  1. 

Clearly,  it  is  important  to  know  which  impurities  are  effective 
at  enhancing  which  native  defect(s).  In  this  article,  we  report  our 
first-principles  studies  of  various  monovalent  (Na,  K,  Cu,  and  Ag), 
divalent  (Mg  and  Zn),  trivalent  (Al),  tetravalent  (Zr,  C,  and  Si),  and 
pentavalent  (V  and  Nb)  impurities  in  LiFeP04  and  possible  defect 
complexes  between  the  extrinsic  impurities  and  native  defects.  The 
effects  of  the  impurities  on  the  ionic  and  electronic  conduction  will 
be  discussed  and  comparisons  with  available  experimental  data 
will  be  made  where  appropriate.  On  the  basis  of  our  studies,  we 
discuss  specific  strategies  for  enhancing  the  electrical  conductivity 
in  LiFeP04. 

2.  Methodology 

Our  first-principles  calculations  of  native  defects  and  extrin¬ 
sic  impurities  (hereafter  referred  to  as  “defects”)  were  based  on 
density-functional  theory  within  the  GGA  +  U  framework  [19-21  ], 
which  is  an  extension  of  the  generalized-gradient  approximation 
[22],  and  the  projector  augmented  wave  method  [23,24],  as  imple¬ 
mented  in  the  VASP  code  [25-27].  The  U  value  for  Fe  was  taken 
from  Zhou  et  al.  [28].  We  used  an  orthorhombic  (1x2x2)  super¬ 
cell,  corresponding  to  112  atoms  per  cell,  to  model  the  defects  in 
LiFeP04. 

Different  defects  in  LiFeP04  are  characterized  by  their  formation 
energies.  The  formation  energy  (£-0  of  a  defect  is  a  crucial  factor  in 
determining  its  concentration.  In  thermal  equilibrium,  the  concen¬ 
tration  of  the  defect  X  at  temperature  T  can  be  obtained  via  the 
relation  [7,18] 

^(X)  =  NsitesNconqg  exp[-£/(X)/kBT],  (1) 

where  Nsites  is  the  number  of  high-symmetry  sites  in  the  lattice  per 
unit  volume  on  which  the  defect  can  be  incorporated,  and  Nconfig 
is  the  number  of  equivalent  configurations  (per  site).  It  emerges 
from  Eq.  (1)  that  defects  with  lower  formation  energies  will  occur 
in  higher  concentrations. 

The  formation  energy  of  a  defect  X  in  charge  state  q  is  defined 
as  [7,18] 

Ef  (XI)  =  Etot(X")  -  Etot(bulk)  -  +  AV  +  /Me),  (2) 

i 

where  Etot(Xq)  and  Etot(bulk)  are  the  total  energies  of  a  supercell 
containing  X  and  of  a  supercell  of  the  perfect  bulk  material;  is 
the  atomic  chemical  potential  of  species  i  (and  is  referenced  to  the 
standard  state),  and  n2  denotes  the  number  of  atoms  of  species  i 
that  have  been  added  (nz> 0)  or  removed  (n,-  <  0)  to  form  the  defect. 
lie  is  the  electron  chemical  potential  (the  Fermi  level)  referenced  to 
the  VBM  in  the  bulk  (Ev).  A  V  is  the  “potential  alignment”  term,  i.e., 
the  shift  in  the  band  positions  due  to  the  presence  of  the  charged 
defect  and  the  neutralizing  background,  obtained  by  aligning  the 
average  electrostatic  potential  in  regions  far  away  from  the  defect 
to  the  bulk  value  [18]. 

The  atomic  chemical  potentials  /x2  are  variables  and  can  be 
chosen  to  represent  experimental  conditions,  and  are  subject  to 
various  thermodynamic  constraints  [7].  These  constraints  ensure 
that  LiFeP04  is  thermodynamically  stable.  The  calculated  phase  dia¬ 
grams  of  the  quaternary  Li— Fe— P— 02  system  at  0  K  reported  by  Ong 
et  al.  [29]  show  that  the  compound  is  stable  over  a  range  of  the 
oxygen  chemical  potential  (/x02)  values,  from  -3.03  to  -8.25  eV. 
These  two  values  are  the  upper  and  lower  limits  in  the  range  of 
fio2  values  considered  in  our  work.  Lower  /x02  values  represent  the 
so-called  “more  reducing  environments,”  which  are  usually  associ¬ 
ated  with  higher  temperatures  and  lower  oxygen  partial  pressures 
and/or  the  presence  of  oxygen  reducing  agents ;  whereas  higher  /xq2 
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values  represent  “less  reducing  environments”  [29].  For  each  /z02, 
the  chemical  potential  of  Li,  Fe,  and  P  can  be  defined  through  a  set 
of  constraints  relevant  at  that  /z02  value.  For  example,  at  /z02  = 
-4.59  eV,  the  thermodynamically  allowed  range  of  /zLi  and  /zFe 
values  is  the  area  defined  by  secondary  phases  Fe203,  Fe3(P04)2, 
Fe2P207,  Li4P207,  and  Li3P04  in  the  chemical-potential  diagram  as 
depicted  in  Fig.  1  of  Ref.  [7];  the  remaining  variable  /zP  is  deter¬ 
mined  via  the  condition  that  ensures  the  stability  of  LiFeP04  [7]. 
Our  investigation  of  thermodynamically  allowed  /zLi,  /z Fe,  /zP,  and 
lio2  values  thus  involves  all  possible  Li— Fe— P— 02  phases  reported 
in  Ref.  [29].  Environments  under  which  the  system  is  close  to  form¬ 
ing  Li-containing  (Fe-containing)  secondary  phases  are  referred  to 
as  Li-excess  (Li-deficient)  [7,29]. 

For  the  extrinsic  impurities,  the  lower  limit  on  /z,  is  minus 
infinity  and  the  upper  limit  is  zero  (with  respect  to  the  energy 
of  the  elemental  bulk  phase)  [18].  Although  stronger  bounds  on 
the  impurity  chemical  potentials  can  be  estimated  based  on  other 
solubility-limiting  phases  formed  between  the  impurities  and  the 
host  constituents,  there  is  no  way  to  know  their  exact  values.  In  the 
following  presentation,  the  atomic  chemical  potentials  of  the  impu¬ 
rities  were  chosen  somewhat  arbitrarily  within  the  range  given  by 
the  above  mentioned  upper  and  lower  limits,  provided  that  the 
calculated  formation  energies  at  /zj,nt  are  positive.  These  choices, 
however,  in  no  way  affect  the  physics  of  what  we  present  since 
we  are  interested  only  in  the  relative  formation  energies  of  the 
impurities  at  a  given  set  of  Li,  Fe,  P,  and  02  chemical  potentials.  In 
other  words,  we  are  trying  to  answer  the  following  question:  Once 
a  certain  impurity  is  incorporated  into  the  material,  what  is  its  lat¬ 
tice  site  preference  under  thermodynamic  equilibrium  conditions 
and  what  is  its  possible  effect  on  the  material’s  properties?  Further 
details  of  the  calculations,  our  methodology  (and  its  limitations), 
and  the  results  for  native  defects  in  LiFeP04  can  be  found  in  Ref. 
[7]. 


3.  Results  and  discussion 

Impurities  in  LiFeP04  can  be  substitutional  at  the  Li,  Fe,  and 
P  sites,  or  may  be  at  interstitial  sites.  For  each  substitutional  or 
interstitial  impurity,  calculations  were  carried  out  in  several  pos¬ 
sible  charge  states.  We  paid  attention  to  positively  and  negatively 
charged  impurities  that  have  the  lowest  formation  energy  at  /Zgnt, 
and  hence  the  highest  concentration,  since  they  can  act  as  acceptor¬ 
like  (donor-like)  dopants  that  shift  the  Fermi  level  toward  the 
VBM  (CBM).  Note  that  /z^nt  is  dependent  on  the  chemical  poten¬ 
tials.  For  example,  we  find  that  /zj,nt  =  0.58  -  0.59  eV  for  /z02  = 
-3.03  eV,  0.72-1.06 eV  for  /z o2  =  -4.59 eV,  and  1.98-2.00 eV  for 
/zo2  =  -8.21  eV  [7].  We  also  considered  possible  complexes  con¬ 
sisting  of  extrinsic  impurities  and  native  point  defects.  In  the 
following,  we  present  only  results  for  those  defect  configurations 
that  have  low  formation  energies  and  hence  are  potentially  relevant 
to  the  material’s  properties.  Note  that  the  defect  formation  ener¬ 
gies  presented  in  all  the  figures  were  evaluated  at  /z02  =  -3.03  eV 
(corresponding  to,  e.g.,  a  temperature  of  800  °C  and  an  oxygen  par¬ 
tial  pressure  of  10-3  atm  [30]),  and  equilibrium  between  LiFeP04, 
Fe203,  and  Fe7(P04)6  was  assumed  (i.e.,  Li-deficient  environment) 
[7].  These  conditions  give  rise  to  /zLi  =-3.41  eV,  /zFe  =  -  3.35  eV,  and 
/zP  =  -  6.03  eV  with  the  formation  enthalpies  of  the  Li— Fe— P— 02 
phases  taken  from  Ref.  [29].  Under  these  conditions,  p+  and  V~ 
are  the  dominant  native  point  defects  and  the  Fermi  level  is  at 
/Zgnt  =  0.59eV  [7],  as  illustrated  in  Fig.  1.  Of  course  the  energy 
landscape  may  be  different  for  a  different  set  of  Li,  Fe,  P,  and  02 
chemical  potentials.  Therefore,  we  will  extend  our  discussions  and 
conclusions  to  all  thermodynamically  allowed  chemical  potentials 
mentioned  in  the  previous  section. 


VBM  CBM 


Fig.  2.  Formation  energies  of  (a)  Na-  and  (b)  K-related  impurities  at  \i  o2  =  -3.03  eV, 
plotted  as  a  function  of  Fermi  level  with  respect  to  the  VBM. 


3.1.  Monovalent  impurities 

Fig.  2(a)  and  (b)  shows  the  formation  energies  of  Na-  and  K- 
related  substitutional  impurities  at  the  Li  and  Fe  sites,  evaluated 
at  /zo2  =  -3.03  eV  and  its  associated  /zLi,  /zFe,  and  /zP  values  as 
noted  earlier.  The  slope  in  the  formation  energy  plots  indicates  the 
charge  state.  Positive  slope  indicates  that  the  defect  is  positively 
charged,  negative  slope  indicates  the  defect  is  negatively  charged; 
cf.  Eq.  (2).  We  find  that  Na°  and  Na“e  are  elementary  point  defects; 
whereas  Naj^  is,  in  fact,  a  complex  of  Na^  and  small  hole  polaron 
p+  that  is  stabilized  at  the  neighboring  Fe  site,  and  Na°e  a  complex 
of  Na“e  and  p+.  At  the  Fe  site  where  the  hole  polaron  is  located, 
the  average  Fe— O  bond  length  is  2.06  A  (compared  to  2.18  A  of  the 
other  Fe— O  bonds)  and  the  calculated  magnetic  moment  is  4.28/zb 
(compared  to  3.76/zB  at  other  Fe  sites);  these  numbers  are  in  excel¬ 
lent  agreement  with  those  obtained  for  an  isolated  hole  polaron  in 
LiFeP04  [7].  For  all  possible  Li,  Fe,  P,  and  02  chemical  potentials, 
Na°  has  the  lowest  formation  energy  at  /zj,nt,  as  seen  in  Fig.  2(a) 
for  /zq2  =  -3.03  eV,  suggesting  that  Na  is  likely  to  exist  as  neutral 
Na°j  under  thermodynamic  equilibrium  and  therefore  not  effec¬ 
tive  in  enhancing  either  p+  or  V~.  The  formation  energy  difference 
between  Na°  and  NajTe  is,  however,  small  (about  0.25  eV  at  most, 
as  measured  at  /z^nt). 

Strictly  speaking,  Na+  is  not  the  positive  charge  state  of  substi¬ 
tutional  Na  at  the  Li  site,  and  Na°e  is  not  the  neutral  charge  state 
of  Na  at  the  Fe  site.  In  other  words,  NaLi,  as  a  single  point  defect,  is 
only  stable  as  Na°,  and  NaFe  as  Naj:e,  since  locally  stable  configura¬ 
tions  of  charge  states  other  than  Na°  and  Na“e  cannot  be  stabilized. 
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If  we  try,  e.g.,  to  create  Na+,  it  decays  to  a  situation  where  the 
positive  charge  is  not  associated  with  the  point  defect  but  corre¬ 
sponds  to  a  small  hole  polaron  stabilized  at  the  neighboring  Fe  site, 
as  presented  above.  Note  that,  instead  of  writing  down,  e.g.,  Na+  as 
Na°-p+,  we  keep  the  notation  Na+  for  simplicity,  but  it  should  be 
regarded  as  a  nominal  notation  referring  to  a  defect  complex  that 
consists  of  the  stable/elementary  Na°  and  a  small  hole  polaron. 

For  K-related  impurities,  I<°  and  K“e  are  elementary  defects; 
whereas  K+  is  a  complex  of  I<°  and  p+,  and  Kpe  a  complex  of  Kpe  and 
p+,  all  similar  to  Na-related  impurities.  Flowever,  we  find  that,  for 
Ho2  =  -3.25  eV  and  Li-excess  or  for  -3.25  eV  <  iiq2  <  -3.03  eV, 
K“e  and  K+  have  the  lowest  energies  for  the  entire  range  of  the 
Fermi-level  values,  similar  to  the  energy  landscape  presented  in 
Fig.  2(b);  whereas  K°  (or  Kpe,  depending  on  the  specific  set  of  the 
chemical  potentials)  has  the  lowest  energy  for  /x02  =  -3.25  eV  and 
Li-deficiency  or  for  -8.21  eV  <  /x02  <  -3.25  eV.  The  fact  that  Kpe 
can  have  lower  formation  energy  than  any  other  K-related  impurity 
at  /4nt  suggests  that  this  defect  can  shift  the  Fermi  level  toward  the 
VBM  and  lower  the  formation  energy  of  positively  charged  native 
defects.  As  a  result,  there  will  be  an  increase  in  the  concentration 
ofp+  and  hence  in  the  electronic  conductivity.  It  is,  however,  noted 
that  the  formation  energy  difference  between  K-  and  K°  (or  Kpe)  is 
small  (about  0.15  eV  at  most).  The  difference  between  the  general 
results  for  Na  and  those  for  K,  albeit  small,  can  be  attributed  to  the 
difference  in  their  ionic  radii.  Note  that  K  can  be  incorporated  at 
the  Fe  site  only  for  a  small  range  of  the  oxygen  chemical  potential 
(-3.25  eV  <  fi o2  <  -3.03  eV)  under  thermodynamic  equilibrium, 
and  the  Fermi-level  shift  is  relatively  small  (about  0.18  eV  at  most). 
The  effects  of  K,  therefore,  may  not  be  very  significant. 

Similar  calculations  were  also  carried  out  for  Cu-  and  Ag-related 
impurities.  We  find  that  Cu°  has  the  lowest  formation  energy  for 
/jLo2  <  -4.59  eV,  whereas  CUpe  has  the  lowest  energy  for  /x02  > 
-4.59  eV.  For  Ag-related  impurities,  Ag°  has  the  lowest  forma¬ 
tion  energy  at  /4nt  f°r  thermodynamically  allowed  Li,  Fe,  P,  and 
02  chemical  potentials.  Therefore,  Cu  and  Ag  are  not  likely  to  be 
effective  in  shifting  the  Fermi  level,  similar  to  Na. 

Experimentally,  it  has  been  suggested  that  Na-  and  K-doped 
LiFeP04/C  samples  have  an  improved  electrochemical  perfor¬ 
mance  compared  to  the  undoped  LiFeP04/C  [31,32].  For  example, 
Na  doping  was  reported  to  slightly  enhance  the  electrical  con¬ 
ductivity;  1.9  x  10-2  and  0.55  x  10-2  Son-1  for,  respectively, 
Li0.97Na0.03FePO4/C  and  undoped  LiFeP04/C  at  room  temperature 
[31].  As  discussed  above,  Na  at  the  Li  site  is  electrically  inactive 
and  thus  cannot  enhance  the  concentration  of  p+  or  V~  and  hence 
the  conductivity.  Furthermore,  Na  is  less  likely  to  be  incorporated 
at  the  Fe  site.  One,  however,  cannot  exclude  the  possibility  of  also 
having  a  significant  concentration  of  Na  at  the  Fe  sites,  especially 
if  the  samples  are  prepared  under  conditions  far  from  equilibrium. 
Given  the  relatively  small  formation  energy  difference  between  the 
Li  and  Fe  sites,  further  experimental  studies  should  be  carried  out 
to  clarify  the  lattice  site  preference  of  the  monovalent  impurities 
in  LiFeP04  and  its  dependence  on  the  synthesis  conditions. 


3.2.  Divalent  impurities 

Fig.  3  shows  the  formation  energies  of  Mg-related  substitutional 
impurities  at  the  Li  and  Fe  sites,  evaluated  at  /x02  =  -3.03  eV.  We 
find  that  Mg°e  and  Mg+  are  elementary  point  defects,  whereas 
Mg^+  is  a  complex  of  Mg+  and  p+.  Among  all  possible  Mg-related 
impurities,  Mg°e  has  the  lowest  formation  energy  at  /4nt  for  all  pos¬ 
sible  Li,  Fe,  P,  and  02  chemical  potentials.  Note  that  the  formation 
energy  at  /4nt  °f  Mg+  *s  higher  than  that  of  Mg°e  by  0.13-0.36  eV, 
depending  on  the  specific  set  of  the  chemical  potentials.  In  fact, 
the  energy  of  Mg°e  is  lowest  for  almost  the  entire  range  of  the 
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Fig.  3.  Formation  energies  of  Mg-related  impurities  at  /x02  =  -3.03  eV,  plotted  as  a 
function  of  Fermi  level  with  respect  to  the  VBM. 

Fermi-level  values  except  very  near  the  VBM,  as  seen  in  Fig.  3  for 
fio2  =  -3.03  eV.  Thus  Mg  is  likely  to  be  electrically  inactive.  Zn 
gives  similar  results  (not  shown  in  the  figure),  i.e.,  Zn  is  likely  to 
exist  as  neutral  Zn°e  and  does  not  shift  /xe  away  from  /4nt. 

Experimentally,  Roberts  et  al.  [33]  reported  that  there  was  no 
evidence  of  Mg  on  the  Li  site  in  samples  prepared  with  the  nom¬ 
inal  stoichiometry  Li!_xMg0.5XFePO4,  which  is  consistent  with  our 
conclusion  that  Mg  is  likely  to  exist  as  Mg°e  in  LiFeP04.  Note  that 
these  results  are  in  contrast  to  those  reported  by  other  research 
groups  who  suggested  that  Mg  occurs  on  the  Li  site  [4,34,35]. 
Since  Mg°e  can  enhance  neither  the  concentration  of  p+  nor  that 
of  V~ ,  the  reported  enhancement  in  the  electrical  conductivity  of 
Mg-doped  LiFeP04  [4,33,36]  is  difficult  to  understand  as  arising 
directly  from  Mg  doping  at  the  Fe  site.  Flowever,  given  the  rela¬ 
tively  small  formation  energy  difference  at  /4nt  between  the  Li  and 
Fe  sites  (0.13-0.36  eV),  a  significant  concentration  of  Mg  at  the  Li 
sites  might  still  be  possible,  especially  if  the  samples  are  prepared 
under  conditions  far  from  equilibrium.  Regarding  Zn,  Bilecka  et  al. 
[37]  reported  a  conductivity  of  10.42xl0-8 Son-1  at  room  tem¬ 
perature  in  LiZnxFei_xP04  samples,  compared  to  7.5 xlO-8  Son-1 
of  undoped  LiFeP04.  Since  Zn°e  is  not  effective  in  shifting  the  Fermi 
level,  this  small  conductivity  increase  may  be  due  to  other  effects. 

3.3.  Trivalent  impurities 

Fig.  4  shows  the  formation  energies  of  Al-related  substitutional 
impurities  at  the  Fe  and  P  sites,  evaluated  at  /x02  =  -3.03  eV.  We 
find  that  Al+e  is  an  elementary  point  defect,  whereas  Alp+  is  a  com¬ 
plex  of  Alpe  and  p+.  At  the  P  site,  Alp-  is  found  to  have  the  lowest 
energy,  which  is  not  unexpected,  given  the  valence  of  Al  (+3)  and 
P  (+5).  In  this  configuration,  Al  replaces  P  and  forms  a  slightly 
distorted  A104  unit  with  the  average  Al— O  bond  length  of  1.77  A, 
compared  to  1.56  A  of  the  P—O  bonds.  Among  the  Al-related  impu¬ 
rities,  we  find  that  A1+  has  the  lowest  formation  energy  at  /4nt  for 
all  thermodynamically  allowed  Li,  Fe,  P,  02  chemical  potentials.  The 
formation  energy  of  the  lowest  energy  configuration  of  AlLi  is  found 
to  be  higher  than  that  of  AlFe  by  0.65-0.89  eV,  as  measured  at  /4nt. 

With  Al  existing  as  the  positively  charged  impurity  Alpe,  more 
V~  and  Li-e  are  created  to  maintain  charge  neutrality,  thus  the 
Fermi  level  is  shifted  away  from  /4nt  and  toward  the  CBM.  As 
a  result  of  this  shift,  the  formation  energy  (concentration)  of 
p+  is  increased  (decreased).  We  also  investigated  possible  defect 
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Fig.  4.  Formation  energies  of  Al-related  impurities  at  /xo2  =  -3.03  eV,  plotted  as  a 
function  of  Fermi  level  with  respect  to  the  VBM. 

complexes  between  Alpe  and  the  negatively  charged  native  point 
defects,  i.e.,Al+e  -  V~  and  Al+-Li“e.  We  find  that  the  binding  energy 
of  these  complexes  are,  respectively,  0.38  and  0.25  eV  with  respect 
to  their  isolated  constituents.  For  any  given  set  of  the  atomic 
chemical  potentials,  Al+-Li“e  has  a  slightly  higher  (0.04-0.23  eV) 
formation  energy  than  Al+-17“.  The  distance  between  the  two 
defects  in  Al+e-V“  is  3.28  A,  whereas  it  is  3.96  A  in  Al+e-Li“e.  Our 
results  suggest  that,  although  Al+e  enhances  both  Vjr  and  Li“e,  the 
concentration  enhancement  is  larger  for  the  vacancies  than  for  the 
antisites. 

Experimentally,  Amin  et  al.  [38,39]  reported  that  Al-doped 
LiFeP04  has  a  higher  ionic  conductivity  but  lower  electronic  con¬ 
ductivity  compared  to  undoped  LiFeP04.  This  is  exactly  consistent 
with  our  results  that  Alpe  enhances  V~  and  reduces  p+.  Regarding 
the  lattice  site  preference,  although  Amin  et  al.  suggested  that  Al 
occurs  on  the  Fe  site,  which  is  in  agreement  with  our  results,  other 
research  groups  reported  that  the  dopant  is  predominantly  on  the 
Li  site  [4,6,35].  Further  experimental  studies,  therefore,  should  be 
carried  out  to  clarify  this  situation. 

3.4.  Tetravalent  impurities 

Fig.  5  shows  the  formation  energies  of  Zr- related  substitutional 
impurities  at  the  Fe  and  P  sites,  evaluated  at /xo2  =  -3.03  eV.  Poten¬ 
tially  relevant  defect  configurations  are  Zr3+,  Zr3+  (a  complex  of 
Zrp+  and  p+),  and  Zr“.  In  Zr“,  Zr  replaces  P  and  forms  a  Zr04  unit 
with  the  average  Zr— O  bond  length  of  1 .98  A,  compared  to  1 .56  A  of 
the  P—O  bonds.  We  find  that  Zrj:+  has  the  lowest  formation  energy 
at  /4nt  for  all  possible  Li,  Fe,  P,  02  chemical  potentials,  suggesting 
that  Zr  has  a  donor-like  effect,  i.e.,  shifting  the  Fermi  level  toward 
CBM  and  hence  enhancing  the  negatively  charged  native  defects. 
The  formation  energy  of  the  lowest  energy  configuration  of  ZrLi  at 
/4nt  is  found  to  be  higher  than  that  of  ZrFe  by  1.17-1.41  eV. 

We  also  considered  complexes  between  the  elementary  Zr2+ 
and  the  negatively  charged  native  point  defects.  Zr2+ -217“  is  found 
to  have  a  binding  energy  of  0.78  eV  with  respect  to  its  isolated  con¬ 
stituents;  the  distances  between  Zr2+  and  the  vacancies  are  3.26 
and  3.62  A.  Zr2+-2Li“e  has  a  binding  energy  of  0.65  eV;  the  dis¬ 
tances  between  Zrj:+  and  the  antisites  are  4.05  and  4.13  A.  Under 
Li-deficient  environments,  the  formation  energy  of  Zr{:+ -217“  is 
found  to  be  lower  than  that  of  Zrj:+ -2Li“e  by  about  0.20-0.34  eV  as 


VBM  CBM 


Fig.  5.  Formation  energies  ofZr-related  impurities  at  jll o2  =  -3.03  eV,  plotted  as  a 
function  of  Fermi  level  with  respect  to  the  VBM. 


seen  in  Fig.  5;  whereas  under  Li-excess  environments,  the  energy  of 
Zr2+ -217“  is  comparable  to  that  of  Zr2+ -2Li“e.  Our  results  therefore 
suggest  that  Zr2+  is  more  likely  to  be  associated  with  17“  than  with 
Li“e,  thus  enhancing  the  vacancies  and  the  ionic  conductivity. 

The  donor-like  effect  of  Zr  in  LiFeP04  could  have  contributed  to 
the  enhancement  in  the  total  conductivity  as  observed  in  exper¬ 
iment  [4].  Flowever,  in  order  to  clearly  see  the  effect,  one  should 
carefully  look  at  the  ionic  and  electronic  conductivities  in  undoped 
and  doped  samples  (preferably  without  carbon  coating).  Note  that, 
an  enhancement  in  the  intrinsic  ionic  conductivity  usually  results 
in  a  reduction  in  the  intrinsic  electronic  conductivity,  as  discussed 
in  the  Introduction.  Regarding  the  lattice  site  preference,  several 
experimental  works  suggested  that  Zr  occupies  the  Li  site  [4,6,35], 
which  is  in  contrast  to  our  results  showing  that  the  Fe  site  is 
energetically  more  favorable  under  thermodynamic  equilibrium 
conditions. 

Fig.  6(a)  and  (b)  shows  the  formation  energies  of  C-  and  Si- 
related  substitutional  impurities  at  the  P  site,  evaluated  at  /x02  = 
-3.03  eV.  Cp  is  an  elementary  defect,  whereas  Cf]  is  a  complex  of 
Cp  and  p+,  Cp  a  complex  of  Cp  and  two  p+.  We  also  considered 
Cp  -17q+,  a  positively  charged  complex  consisting  of  Cp  and  an  oxy¬ 
gen  vacancy  (I7q+).  This  complex  has  a  binding  energy  of  4.44  eV 
with  respect  to  its  isolated  constituents.  It  can  also  be  regarded  as 
replacing  a  (P04)3-  tetrahedron  with  a  trigonal  planar  (C03)2-.  In 
the  (C03  )2“  unit,  the  average  C— O  bond  length  is  1.31  A,  compared 
to  1 .56  A  of  the  P—O  bonds.  C  thus  exhibits  a  +4  valence  at  the  P  site. 
At  the  Li  and  Fe  sites,  we  however  find  that  C  exhibits  a  +2  valence. 
As  a  result,  C+  and  C|?e  are  the  elementary  defects,  whereas  their 
higher  charge  states  are  defect  complexes  consisting  of  C+  (or  C°e) 
and  p+. 

For  -3.25  eV  <  /x02  <  -3.03  eV,  Cjt  has  the  lowest  formation 
energy  at  /4nt.  Given  its  positive  effective  charge,  this  defect  would 
exhibit  a  donor-like  effect,  enhancing  the  concentration  of  nega¬ 
tively  charged  native  defects.  Moreover,  for  each  C^  created,  there 
are  two  p+;  thus,  with  Cp  being  the  lowest  energy  configuration, 
C  can  enhance  both  p+  and  V~  (see  more  below).  For  -8.21  eV  < 
Ho2  <  -3.25  eV,  the  formation  energy  of  C“  -  172+  is  lower  than 
that  of  C+,  and  the  complex  has  the  lowest  energy  at  /4nt.  Since 
Cp -17q+  exhibits  a  donor-like  effect,  it  can  enhance  V~ .  In  fact,  we 
have  investigated  a  complex  of  Cp  -172+  and  V~ ,  hereafter  denoted 
as  Cp  -  17q+  -  \7“,  and  find  that  it  has  a  binding  energy  of  0.63  eV 
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Fig.  6.  Formation  energies  of  (a)  C-  and  (b)  Si-related  impurities  at  /jl02  =  -3.03  eV, 
plotted  as  a  function  of  Fermi  level  with  respect  to  the  VBM. 

with  respect  to  its  isolated  constituents.  The  formation  energy  of 
this  complex  [see  Fig.  6(a)]  is  lower  than  that  of  Cp -V“  (not  shown 
in  the  figure),  a  neutral  complex  of  and  V~,  by  0.25-3.79  eV. 
Our  results  suggest  that,  in  the  presence  of  V£7,  C“  -V2+  is  ener¬ 
getically  more  favorable  than  Cp  for  all  possible  Li,  Fe,  P,  and  02 
chemical  potentials.  In  other  words,  if  incorporated  into  the  bulk 
of  LiFeP04,  C  is  likely  to  exist  as  Cp  and  the  presence  of  C 

thus  enhances  V~  only.  The  formation  energy  of  C  at  the  Li  and  Fe 
sites  is  much  higher  than  at  the  P  site,  with  the  energy  difference  at 
/4nt  between  Cp  and  Cpe  is  4.80-6.24  eV.  Note  that,  here  we  assume 
C  can  go  into  the  bulk  with  a  concentration  higher  than  that  of  the 
native  defects.  The  penetration  of  C,  however,  may  be  impeded  by 
reactions  at  the  surface  [40].  Since  C  is  present  in  precursors  used  in 
the  synthesis  of  LiFeP04  and  also  because  LiFeP04  is  usually  coated 
with  C,  it  is  worth  carrying  out  further  studies  to  see  if  there  is  C  in 
the  interior  of  the  material. 

Regarding  Si-related  substitutional  impurities  at  the  P  site,  Sip 
is  found  to  be  an  elementary  defect,  whereas  Sip  is  a  complex  of 
Si“  and  p+,  and  Si+  a  complex  of  Si“  and  two  p+.  Sip-V2+  has 
a  binding  energy  of  2.35  eV  with  respect  to  Sip  and  V^+.  In  this 
complex,  Si  forms  a  Si04  unit  that  shares  one  O  atom  with  the  neigh¬ 
boring  P04  unit,  i.e.,  Si03— O— P03  with  the  Si— P  distance  being 
2.90  A,  compared  to  3.71  A  of  the  P— P  distance.  Energetically,  we 
find  that  Sip  has  a  lower  formation  energy  than  Sip-Vo+  except 
for  -8.21  eV  <  /x02  <  -7.59  eV.  Sij]  is  found  to  have  the  lowest  for¬ 
mation  energy  at  /4nt  f°r  -4.59 eV  <  fi0  <  -3.03  eV,  as  seen  in 
Fig.  6(b)  for  /x02  =  -3.03  eV.  Although  Si"  is  neutral  and  thus  not 
effective  in  shifting  the  Fermi  level,  each  Sip  contains  one  p+.  The 


formation  of  Sip,  therefore,  may  help  introduce  more  p+  into  the 
material.  For  -8.21  eV  <  /x02  <  -4.59  eV,  Sip  is  found  to  have  the 
lowest  formation  energy.  In  this  case,  Si  exhibits  an  acceptor-like 
effect  and  is  expected  to  enhance  the  concentration  of  p+.  For  any 
given  set  of  Li,  Fe,  P,  and  02  chemical  potentials,  neutral  complexes 
such  as  those  consisting  of  V~  and  Si+  or  Si“ -V2+  all  have  higher 
formation  energies  than  Si°.  Our  results  thus  suggest  that  Si  incor¬ 
porated  at  the  P  site  enhances  small  hole  polarons.  At  the  Li  and  Fe 
sites,  Si  exhibits  a  donor-like  effect  since  Si  is  supervalent  to  Li  and 
Fe;  the  elementary,  lowest  energy  configurations  at  /4nt  are  Si3+ 
and  Sip+,  where  the  formation  energy  of  Si  at  the  Li  site  is  higher 
than  that  at  the  Fe  site  by  0.43-0.65  eV,  depending  on  the  specific 
set  of  the  chemical  potentials.  The  formation  energy  difference  at 
/4nt  between  SiP  and  SiFe  is  0.33-3.56  eV 

Amin  et  al.  [41  ]  reported  that  doping  LiFeP04  with  Si  results  in 
an  increase  of  the  ionic  conductivity  and  a  decrease  of  the  elec¬ 
tronic  conductivity,  suggesting  a  donor-like  effect.  As  discussed 
above,  this  effect  should  be  associated  with  Si2+.  Our  results,  how¬ 
ever,  suggest  that  Sij:+  is  not  the  lowest  defect  configuration  under 
thermodynamic  equilibrium  conditions.  Although  the  cause  of  this 
discrepancy  between  experiment  and  our  results  on  the  lattice  site 
preference  of  Si  is  still  not  clear,  we  suspect  that  doping  at  the  P  site 
may  be  kinetically  hindered,  and  thus  thermodynamic  equilibrium 
approaches  have  limited  applicability  [18].  This  may  be  the  case 
for  all  dopants  that  involve  P-site  doping  in  LiFeP04.  It  should  be 
noted  that  LiFeP04  can  be  regarded  as  consisting  of  Li+,  Fe2+,  and 
(P04)3-  with  ionic  bonding  between  the  units.  Within  the  (P04)3- 
unit,  however,  the  P—O  bonds  are  highly  covalent  [7].  If  doping 
occurs  via  diffusion,  P-site  doping  would  involve  diffusion  of  P5+ 
and  the  dopant  which  is  expected  to  be  difficult  because  it  neces¬ 
sarily  involves  breaking/forming  of  chemical  bonds  between  O  and 
P  (and  the  dopant).  Thus  P-site  doping  might  also  be  affected  by  the 
P-  and  dopant-containing  precursors  used  in  the  synthesis. 

3.5.  Pentavalent  impurities 

Fig.  7(a)  and  (b)  shows  the  formation  energies  of  V-  and  Nb- 
related  substitutional  impurities  at  the  Li,  Fe,  and  P  sites,  evaluated 
at  fio2  =  -3.03  eV.  We  find  that  V°  is  an  elementary  point  defect, 
suggesting  that  V  (vanadium,  not  to  be  confused  with  a  vacancy, 
V)  at  the  P  site  has  a  +5  valence.  V+  is  a  complex  of  V°  and  p+, 
whereas  V“  is  a  complex  of  Vp  and  a  small  electron  polaron  (p_) 
[7].  At  the  Li  and  Fe  sites,  however,  V  exhibits  a  +3  valence,  with 
V2+  and  Vpe  being  the  elementary  defects.  Attempts  to  create  higher 
charge  states  result  in  complexes  consisting  of  these  defects  and  p+, 
whereas  lower  charge  states  are  complexes  containing  p~ .  We  have 
also  investigated  neutral  complexes  such  as  V+-Vu,  a  complex  of 
V+e  and  a  Li  vacancy,  and  2V+ -V2e-,  a  complex  of  two  V+e  and  an 
Fe  vacancy.  Energetically,  we  find  that  Vp  has  the  lowest  formation 
energy  at  /4nt  for  all  possible  Li,  Fe,  P,  and  02  chemical  potentials. 
With  this  neutral  configuration,  V  would  be  ineffective  in  shifting 
the  Fermi  level. 

V  has  positive  effective  charges  at  the  Li  and  Fe  sites  and  thus 
exhibits  a  donor-like  effect.  The  formation  energy  at  /4nt  of  VFe 
in  its  lowest  energy  configuration  is  higher  than  that  of  V°  by 
0.50-2.14  eV  for  almost  all  possible  chemical  potentials,  and  by 
0.1 0-0.20  eV  for  a  small  region  of  the  chemical  potentials  associated 
with  /xo2  around  -4.59  eV.  On  the  other  hand,  the  formation  energy 
of  VFe  is  lower  than  that  of  VLi  by  0.76-0.97  eV  for  almost  all  pos¬ 
sible  chemical  potentials,  except  for  a  small  region  associated  with 
lio2  <  -4.59  eV  where  the  two  defects  have  comparable  formation 
energies.  Neutral  complexes  such  as  Vj^-V”  and  2V+e-V2e_  have 
formation  energies  higher  than  V°  by  0.36-2.32  and  4.51 -9.09  eV, 
respectively.  Again,  these  results  suggest  that  V  is  energetically 
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Fig.  7.  Formation  energies  of  (a)  V-  and  (b)  Nb-related  impurities  at  y io2  =  -3.03  eV, 
plotted  as  a  function  of  Fermi  level  with  respect  to  the  VBM.  For  substitutional 
impurities  at  the  Li  and  Fe  sites  (M^+  and  ,  where  M  =  V  and  Nb,  and  q> 0),  only 
the  lowest  formation-energy  segments  are  presented. 

more  favorable  at  the  P  site  under  thermodynamic  equilibrium 
conditions. 

Experimentally,  Hong  et  al.  [42]  reported  that  their  attempt  to 
incorporate  V  into  LiFeP04  at  the  Fe  site  resulted  in  having  V  at 
the  P  site  instead.  However,  Omenya  et  al.  [43]  later  reported  that 
the  substitution  at  the  P  site  could  not  be  reproduced,  but  at  least 
10  mol%  of  the  Fe  sites  were  occupied  by  V3+.  Zhang  et  al.  [44]  also 
reported  a  valence  between  +3  and  +4  for  V  in  LiFeP04.  Note  that, 
according  our  results  discussed  above,  the  +3  valence  of  V  is  associ¬ 
ated  with  substitution  at  the  Fe  site  that  exhibits  a  donor-like  effect. 
Thus,  for  V-doped  LiFeP04  samples  reported  in  Refs.  [43,44],  V  may 
indeed  have  been  incorporated  at  the  Fe  site,  although  this  site  is 
found  to  be  energetically  less  favorable  in  our  calculations.  The  dis¬ 
crepancy  between  experiment  and  our  results  on  the  lattice  site 
preference  of  V  could  be  that  P-site  doping  is  thermodynamically 
favorable  but  kinetically  hindered,  as  discussed  earlier. 

Regarding  Nb,  we  find  that  Nb°  is  an  elementary  point  defect, 
and  the  dopant  thus  has  a  +5  valence  at  the  P  site.  Nb+  is  a  com¬ 
plex  of  Nb£  and  p+,  whereas  Nb“  is  a  complex  of  Nbf]  and  p~ .  Nb 
also  exhibits  the  +5  valence  at  the  Fe  site,  with  Nbp+  being  the  ele¬ 
mentary  defect.  At  the  Li  site,  we  however  find  that  Nb  has  a  +4 
valence  and  the  elementary  configuration  is  Nb3+.  Higher  charge 
states  such  as  Nb^  is  a  complex  of  Nb3+  and  p+.  Energetically, 
Nbp  has  the  lowest  formation  energy  at  /4nt  for  almost  all  possi¬ 
ble  chemical  potentials,  except  for  a  very  small  region  associated 
with  -4.59  eV  <  /xq2  <  -3.89  eV  of  the  (/xLi,  /xFe,  Ho2 )  polyhedron 


that  defines  the  stable  range  of  Li,  Fe,  P,  and  02  chemical  potentials 
[7]  where  Nbj?+  has  the  lowest  formation  energy.  For  example,  in 
the  pio2  =  -4.59eV  plane,  the  small  region  near  Points  B,  C,  and  D 
in  the  chemical-potential  diagram  as  depicted  in  Fig.  1  of  Ref.  [7] 
is  where  Nbp+  is  most  favorable  energetically.  Nbp  cannot  enhance 
the  concentration  of  either  V~  orp+,  but  Nbj?+  is  donor-like  and  thus 
can  enhance  V~.  We  also  find  that  Nbp+-3VIr,  a  complex  of  Nb3+ 
and  three  Li  vacancies,  has  a  lower  formation  energy  than  Nb°  in 
that  small  region  of  the  (/xLi,  /xFe,  /^o2 )  polyhedron.  Our  results  thus 
suggest  that  Nb  can  be  incorporated  at  the  Fe  or  P  site,  depending 
on  the  synthesis  conditions.  The  reported  enhancement  in  the  total 
conductivity  of  Nb-doped  LiFeP04  [4]  could  have  been  partly  due 
to  the  enhancement  of  V~  caused  by  Nb3+.  In  addition,  the  incorpo¬ 
ration  of  Nb  at  the  P  site  may  be  kinetically  hindered  as  discussed 
above  in  the  case  of  Si  and  V. 

Overall,  we  find  that  extrinsic  impurities  in  LiFeP04  each  have 
one  stable  charge  state  (i.e.,  elementary  defect  configuration)  at  a 
given  lattice  site.  Any  attempt  to  create  higher  (lower)  charge  states 
will  result  in  complexes  consisting  of  the  elementary  configuration 
and  small  hole  (electron)  polarons,  which  is  similar  to  what  has 
been  observed  for  native  defects  in  LiFeP04  as  reported  in  Ref.  [7]. 
Our  studies  help  identify  specific  impurities  that  are  likely  to  be 
effective  (or  ineffective)  for  enhancing  the  concentration  of  p+  or 
V~.  The  results  also  suggest  further  experiments  to  investigate  the 
lattice  site  preference  (and  its  dependence  on  the  synthesis  con¬ 
ditions)  and  changes  in  the  ionic  and  electronic  conductivities.  For 
practical  applications,  since  doping  alone  cannot  enhance  both  p+ 
and  V~  as  needed  for  high  ionic  and  electronic  conductivities,  other 
methods  such  as  carbon  coating  [45,46]  and/or  thermal  treatment 
[47]  may  still  be  needed.  As  discussed  by  Julien  et  al.  [40],  carbon 
coating  can  help  repair  structural  damage  at  the  surface  and  ensure 
the  electric  contact  between  LiFeP04  particles.  Also,  with  whatever 
method  is  employed  to  enhance  the  conductivity,  one  should  first 
be  able  to  control  the  experimental  conditions  during  synthesis  to 
reduce  Fe+  and  enhance  p+  and  V~  [7]. 

Note  that  computational  studies  of  impurities  in  LiFeP04  have 
also  been  carried  out  by  Islam  et  al.  [48,49].  Based  on  calcula¬ 
tions  using  interatomic  potentials,  they  reported  that,  among  other 
impurities,  Na  is  energetically  more  favorable  at  the  Li  site,  whereas 
Mg,  Al,  Zr,  and  Nb  are  more  favorable  at  the  Fe  site.  These  results 
are  generally  in  qualitative  agreement  with  our  results  for  these 
substitutional  impurities  at  the  Li  and  Fe  sites.  However,  they  did 
not  report  results  for  the  impurities  at  the  P  site,  and  also  did  not 
explore  all  possible  Li,  Fe,  P,  and  02  chemical  potentials.  As  dis¬ 
cussed  earlier,  the  lattice  site  preference  can  be  sensitive  to  the 
chemical  potentials  (e.g.,  in  the  case  of  K  and  Cu),  and  the  Li  and 
Fe  sites  may  not  be  energetically  most  favorable  under  thermody¬ 
namic  equilibrium. 

4.  Summary 

We  have  carried  out  first-principles  studies  of  the  effects  of 
various  extrinsic  impurities  on  the  ionic  and  electronic  conduc¬ 
tion  in  LiFeP04.  We  find  that  the  formation  energy  and  lattice  site 
preference  of  the  impurities  depend  on  Li,  Fe,  P,  and  02  chemical 
potentials  which  represent  the  synthesis  conditions.  For  all  ther¬ 
modynamically  allowed  atomic  chemical  potentials,  Na,  Cu,  Ag,  Mg, 
and  Zn  are  likely  to  exist  as  neutral  defects  in  LiFeP04  under  ther¬ 
modynamic  equilibrium  and  thus  do  not  enhance  the  concentration 
of  either  small  hole  polarons  (p+)  or  lithium  vacancies  (!/“),  i.e. 
electronic  or  ionic  conduction.  I<  may  be  incorporated  at  the  Fe  site 
for  a  small  range  of  the  chemical  potentials,  where  it  exhibits  an 
acceptor-like  effect  and  thus  can  enhance  the  concentration  ofp+. 
Al  and  Zr,  on  the  other  hand,  exhibit  a  donor-like  effect  and  are  thus 
effective  in  enhancing  the  concentration  of  V~.  C  can  also  exhibit 
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a  donor-like  effect  if  incorporated  into  the  material.  We  find  that 
Si  is  energetically  more  favorable  at  the  P  site  where  it  exhibits  an 
acceptor-like  effect,  which  is  in  contrast  to  experiment  suggesting 
substitution  at  the  Fe  site  that  exhibits  a  donor-like  effect.  V  is  found 
to  be  more  favorable  at  the  P  site  and  thus  electrically  inactive,  in 
contrast  to  experiment  suggesting  substitution  at  the  Fe  site.  Nb 
is  also  found  to  be  electrically  inactive  at  the  P  site,  except  for  a 
small  range  of  the  atomic  chemical  potentials  where  it  is  energeti¬ 
cally  more  favorable  at  the  Fe  site  and  exhibits  a  donor-like  effect. 
We  suggest  that  doping  LiFeP04  with  Si,  V,  or  Nb  at  the  P  site  is 
thermodynamically  favorable  but  kinetically  hindered,  and  that  Si, 
V,  and  Nb  may  have  been  incorporated  at  the  Fe  site  where  they 
exhibit  a  donor-like  effect  as  reported  in  experiment.  Indeed,  the 
valence  reported  for  V  in  V-doped  LiFeP04  by  experiment  is  consis¬ 
tent  with  substitution  at  the  Fe  site  in  our  calculations.  Our  studies, 
therefore,  can  serve  as  guidelines  for  experiment  on  which  lattice 
site  the  impurities  are  most  likely  to  be  incorporated  under  ther¬ 
modynamic  equilibrium  conditions,  and  how  the  impurities  may 
affect  the  ionic  and  electronic  conduction.  To  achieve  a  high  elec¬ 
trical  conductivity,  however,  one  should  combine  the  incorporation 
of  electrically  active  impurities  with  other  methods  such  as  defect- 
controlled  synthesis,  carbon  coating,  and/or  thermal  treatment. 

Acknowledgements 

We  acknowledge  helpful  discussions  with  S.C.  Erwin,  C.S.  Hell- 
berg,  and  J.  Allen,  and  the  use  of  computing  facilities  at  the  DoD  FIPC 
Centers.  K.H.  was  supported  by  the  U.S.  Naval  Research  Laboratory 
through  Grant  No.  NRL-N00173-08-G001,  and  M.D.J.  by  the  Office 
of  Naval  Research. 

References 

[1]  A.K.  Padhi,  K.S.  Nanjundaswamy,  J.B.  Goodenough,  J.  Electrochem.  Soc.  144 
(1997)1188-1194. 

[2]  B.L.  Ellis,  K.T.  Lee,  L.F.  Nazar,  Chem.  Mater.  22  (2010)  691-714. 

[3]  A.  Manthiram,  J.  Phys.  Chem.  Lett.  2  (2011)  176-184. 

[4]  S.  Chung,  J.  Bloking,  Y.  Chiang,  Nat.  Mater.  1  (2002)  123-128. 

[5]  N.  Ravet,  A.  Abouimrane,  M.  Armand,  Nat.  Mater.  2  (2003)  702. 

[6]  M.  Wagemaker,  B.L.  Ellis,  D.  Liitzenkirchen-Hecht,  F.M.  Mulder,  L.F.  Nazar, 
Chem.  Mater.  20  (2008)  6313-6315. 

[7]  K.  Hoang,  M.  Johannes,  Chem.  Mater.  23  (201 1 )  3003-3013. 

[8]  P.  Axmann,  C.  Stinner,  M.  Wohlfahrt-Mehrens,  A.  Mauger,  F.  Gendron,  C.M. 
Julien,  Chem.  Mater.  21  (2009)  1636-1644. 

[9]  S.-Y.  Chung,  Y.-M.  Kim,  S.-Y.  Choi,  Adv.  Funct.  Mater.  20  (2010)  4219-4232. 

[10]  F.  Zhou,  K.  Kang,  T.  Maxisch,  G.  Ceder,  D.  Morgan,  Solid  State  Commun.  132 

(2004)181-186. 


[11]  T.  Maxisch,  F.  Zhou,  G.  Ceder,  Phys.  Rev.  B  73  (2006)  104301. 

[12]  B.  Ellis,  L.K.  Perry,  D.H.  Ryan,  L.F.  Nazar,  J.  Am.  Chem.  Soc.  128  (2006) 
11416-11422. 

[13]  K.  Zaghib,  A.  Mauger,  J.B.  Goodenough,  F.  Gendron,  C.M.  Julien,  Chem.  Mater. 
19(2007)3740-3747. 

[14]  A.  Peles,  C.G.  Van  de  Walle,  Phys.  Rev.  B  76  (2007)  214101. 

[15]  K.  Hoang,  C.G.  Van  de  Walle,  Phys.  Rev.  B  80  (2009)  214109. 

[16]  G.B.  Wilson-Short,  A.  Janotti,  K.  Hoang,  A.  Peles,  C.G.  Van  de  Walle,  Phys.  Rev.  B 
80(2009)224102. 

[17]  K.  Hoang,  A.  Janotti,  C.G.  Van  de  Walle,  Angew.  Chem.  Int.  Ed.  50  (2011) 
10170-10173. 

[18]  C.G.  Van  de  Walle,  J.  Neugebauer,  J.  Appl.  Phys.  95  (2004)  3851-3879. 

[19]  V.I.  Anisimov,  J.  Zaanen,  O.K.  Andersen,  Phys.  Rev.  B  44  (1991)  943-954. 

[20]  V.I.  Anisimov,  I.V.  Solovyev,  M.A.  Korotin,  M.T.  Czyzyk,  G.A.  Sawatzky,  Phys.  Rev. 
B 48 (1993) 16929-16934. 

[21  ]  A.I.  Liechtenstein,  V.I.  Anisimov, J.  Zaanen,  Phys.  Rev.  B  52  (1995)  R5467-R5470. 

[22]  J.P.  Perdew,  K.  Burke,  M.  Ernzerhof,  Phys.  Rev.  Lett.  77  (1996)  3865-3868. 

[23]  P.E.  Blochl,  Phys.  Rev.  B  50  (1994)  17953-17979. 

[24]  G.  Kresse,  D.  Joubert,  Phys.  Rev.  B  59  (1999)  1758-1775. 

[25]  G.  Kresse,  J.  Hafner,  Phys.  Rev.  B  47  (1993)  558-561. 

[26]  G.  Kresse,  J.  Furthmiiller,  Phys.  Rev.  B  54  (1996)  11169-11186. 

[27]  G.  Kresse,  J.  Furthmiiller,  Comput.  Mater.  Sci.  6  (1996)  15-50. 

[28]  F.  Zhou,  M.  Cococcioni,  C.  Marianetti,  D.  Morgan,  G.  Ceder,  Phys.  Rev.  B  70 (2004) 
235121. 

[29]  P.S.  Ong,  L.  Wang,  B.  Kang,  G.  Ceder,  Chem.  Mater.  20  (2008)  1798-1807. 

[30]  K.  Reuter,  M.  Scheffler,  Phys.  Rev.  B  65  (2001)  035406. 

[31]  X.  Yin,  K.  Huang,  S.  Liu,  H.  Wang,  H.  Wang,  J.  Power  Sources  195  (2010) 
4308-4312. 

[32]  X.  Fang,  J.  Li,  K.  Huang,  S.  Liu,  C.  Huang,  S.  Zhuang,  J.  Zhang,  J.  Solid  State  Elec¬ 
trochem.  (2011),  doi:10.1007/sl0008-01 1-1426-4. 

[33]  M.R.  Roberts,  G.  Vitins,  J.R.  Owen,  J.  Power  Sources  179  (2008)  754-762. 

[34]  Z.  Guo,  H.  Liu,  S.  Bewlay,  H.  Liu,  S.  Dou,  Synth.  Met.  153  (2005)  113-116. 

[35]  N.  Meethong,  Y.-H.  Kao,  S.A.  Speakman,  Y.-M.  Chiang,  Adv.  Funct.  Mater.  19 
(2009) 1060-1070. 

[36]  C.  Wang,  J.  Hong,  Electrochem.  Solid-State  Lett.  10  (2007)  A65-A69. 

[37]  I.  Bilecka, A.  Hintennach,  M.D.  Rossell,  D. Xie,  P.  Novak,  M.  Niederberger, J.  Mater. 
Chem.  21  (2011) 5881-5890. 

[38]  R.  Amin,  C.  Lin,  J.  Maier,  Phys.  Chem.  Chem.  Phys.  10  (2008)  3519-3523. 

[39]  R.  Amin,  C.  Lin,  J.  Maier,  Phys.  Chem.  Chem.  Phys.  10  (2008)  3524-3529. 

[40]  C.M.  Julien,  A.  Mauger,  K.  Zaghib,  J.  Mater.  Chem.  21  (2011)  9955-9968. 

[41  ]  R.  Amin,  C.  Lin,  J.  Peng,  K.  Weichert,  T.  Acartrk,  U.  Starke,  J.  Maier,  Adv.  Funct. 
Mater.  19  (2009)  1697-1704. 

[42]  J.  Hong,  C.S.  Wang,  X.  Chen,  S.  Upreti,  M.S.  Whittingham,  Electrochem.  Solid- 
State  Lett.  12  (2009)  A33-A38. 

[43]  F.  Omenya,  N.A.  Chernova,  S.  Upreti,  P.Y.  Zavalij,  K.-W.  Nam,  X.-Q.  Yang,  M.S. 
Whittingham,  Chem.  Mater.  23  (201 1 )  4733-4740. 

[44]  L.-L.  Zhang,  G.  Liang,  A.  Ignatov,  M.C.  Croft,  X.-Q.  Xiong,  I.-M.  Hung,  Y.- 
H.  Huang,  X.-L.  Hu,  W.-X.  Zhang,  Y.-L.  Peng,  J.  Phys.  Chem.  C  115  (2011) 
13520-13527. 

[45]  H.  Huang,  S.-C.  Yin,  L.F.  Nazar,  Electrochem.  Solid-State  Lett.  4  (2001) 
A170-A172. 

[46]  B.  Ellis,  W.H.  Kan,  W.R.M.  Makahnouk,  L.F.  Nazar,  J.  Mater.  Chem.  17  (2007) 
3248-3254. 

[47]  R.  Amin,  J.  Maier,  Solid  State  Ionics  178  (2008)  1831-1836. 

[48]  M.S.  Islam,  D.J.  Driscoll,  C.A.J.  Fisher,  P.R.  Slater,  Chem.  Mater.  17  (2005) 
5085-5092. 

[49]  C.A.J.  Fisher,  V.M.H.  Prieto,  M.S.  Islam,  Chem.  Mater.  20  (2008)  5907-5915. 


